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Abstract: The Pierre Auger Observatory has been collecting data in a stable manner since 
January 2004. We present here a study of the cosmic ray composition using events recorded 
in hybrid mode during the first years of data taking. These are air showers observed by 
the fluorescence detector as well as the surface detector, so the depth of shower maximum, 
A max , is measured directly. The cosmic ray composition is studied in different energy ranges by 
comparing the observed average A max with predictions from air shower simulations for different 
nuclei. The change of (X max ) with energy (elongation rate) is used to derive estimates of the 
change in primary composition. 
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Introduction 

Ultra-high-energy cosmic rays are presumed 
to be of extragalactic origin. With increas- 
ing energies, and thus Larmor radii, the galac- 
tic charged particles can not be confined in 
our Galaxy and the galactic cosmic ray accel- 
erator candidates are expected to reach their 
maximum energy well below 10 18 eV. More- 
over, there are no experimental signs of an 
anisotropy of the cosmic ray arrival direction 
at these energies. 

The transition between galactic and extra- 
galactic cosmic rays is therefore believed to 
happen between 10 18 and 10 19 eV where a spec- 
tral break in the cosmic ray flux known as the 
'ankle' or 'dip' is observed. The exact position 
and nature of the transition is still disputed 
and it seems clear that a combined precise mea- 
surement of the particle flux and composition 
in this energy range is needed to be able to 
distinguish between different models of the ex- 
tragalactic cosmic ray component (see [1] for 
recent discussions on this subject). 
For fluorescence detectors (FDs), the observ- 
able most sensitive to the composition is the 
slant depth position A max at which the maxi- 



mum of the longitudinal shower profile occurs. 
Its average value (A max ) at a certain energy E 
is related to the mean logarithmic mass (In A) 
via 

(X max ) = D p [In (E/Eo) - (In A}} + c p , (1) 

where D p denotes the 'elongation rate' [2] of a 
proton, and c p is the average depth of a proton 
with reference energy Eq. Both, D p and c p , 
depend on the nature of hadronic interactions. 
The width of the A max distribution is another 
composition sensitive parameter, since heavy 
nuclei are expected to produce smaller shower- 
to-shower fluctuations than protons. 



Data Analysis 

In this analysis we use hybrid events collected 
by the Pierre Auger Observatory between the 
1st of December 2004 and the 30th of April 
2007. These are showers observed by at least 
one FD and with at least one triggered tank 
recorded by the surface detector. 
In order to ensure a good A max resolution at 
the 20 gem -2 level [3], the following quality 
cuts were applied to the data: The recon- 
structed X max should lie within the observed 
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Figure 1 : Illustration of the effect of the field 
of view of the fluorescence detector on the se- 
lected A" max distribution. Filled areas indicate 
slant depths, which are de-selected by the qual- 
ity cuts. 

shower profile and the reduced \ 2 of a fit with 
a Gaisser-Hillas function [4] should not exceed 
2.5. Moreover, insignificant shower maxima 
are rejected by requiring that the \ 2 °f a nn_ 
ear fit to the longitudinal profile exceeds the 
Gaisser-Hillas fit x 2 by at least four. Finally, 
the estimated uncertainties of the shower max- 
imum and total energy must be smaller than 
40 gem -2 and 20%, respectively. 

In addition, a set of fiducial volume cuts is ap- 
plied to allow for an unbiased measurement of 
the X max -distribution: Energy dependent cuts 
on the zenith angle and the maximum tank- 
core distance ensure a single-tank trigger prob- 
ability near one for protons and iron at all en- 
ergies. 

In order to minimise systematic uncertain- 
ties from the relative timing between the flu- 
orescence and surface detectors, the minimum 
viewing angle under which a shower was ob- 
served is required to be larger than 20°. This 
cut also removes events with a large fraction of 
direct Cherenkov light. 

Moreover, a minimisation of the effect of the 
field of view boundaries of the FDs is of utmost 
importance: The current fluorescence detec- 
tors cover an elevation range from f^i = 1.5° to 
Q.2 = 30° and therefore the observable heights 
for vertical tracks are between R tan Qi < h v < 
R tan H2 , where R denotes the distance of the 
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Figure 2: Dependence of the average mea- 
sured AT max on the upper viewable slant depth 
boundary for showers with energies between 
10 18 and 10 18 - 25 eV. The arrow indicates the 
cut corresponding to an estimated contained 
event fraction of > 95%. 

shower core to the fluorescence detector. That 
is, the farther away from a fluorescence detec- 
tor a track is detected, the smaller becomes the 
observable upper slant depth boundary X up . 
Similarly the lower slant depth boundary X\ ow 
becomes larger for near showers. 
Since in the quality selection it is required 
that the X max is detected within the field of 
view, these slant depth boundaries can severely 
bias the selected X max -distributions, as it is 
sketched in Fig. 1. This bias can be avoided 
by selecting only tracks with geometries corre- 
sponding to an X up -X\ ow range, which is large 
enough to contain most of the parent X max - 
distribution. Therefore, we investigate the de- 
pendence of (X max ) on the field of view bound- 
aries and place fiducial volume cuts at the X up 
and X\ ow values, where the (X max ) starts to 
be constant. An example of this procedure is 
shown in Fig. 2. 

Systematic Uncertainties 

The effect of atmospheric uncertainties on 
the measurement of the shower maximum 
is discussed in detail in [5]. The dominat- 
ing contribution is the long-term validity 
of the monthly average molecular profiles 
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Figure 3: (A" max ) as a function of energy compared to predictions from hadronic interaction models. 
The dashed line denotes a fit with two constant elongation rates and a break-point. Event numbers 
are indicated below each data point. 



used in this analysis, which we estimate 
to be < 6 gem" 2 . Using a full detector 
and atmosphere simulation [6], the profile 
reconstruction algorithm [7] was found to be 
unbiased within 5 gem -2 at all energies. The 
effect of multiple-scattered fluorescence and 
Cherenkov light was estimated to contribute 
about 5 gem -2 by comparing different light 
collection algorithms. 

Re-reconstructing showers with the geometry 
determined from the surface detector data 
alone yields an upper bound on the geometri- 
cal uncertainty of < 6 gem . 
The geometrical bias due to the camera 
alignment uncertainty is below 3 gem -2 and 
the residual acceptance difference [8] between 
proton and iron showers contributes around 
10 gcm~ 2 at lowest energies vanishing rapidly 
to zero above 10 18 eV. 

The total uncertainty is thus around 
< 15 gcm~ 2 at low energies and < 11 gcm~ 2 
above 10 18 eV. Note that in addition the 



current uncertainty of the FD energy scale of 
22% [3] needs to be taken into account. 

Results 

After all cuts are applied, 4329 events remain 
for the composition analysis. In Fig. 3 the 
mean A max as a function of energy is shown 
along with predictions from air shower simula- 
tions [10,11]. As can be seen, our measurement 
favours a mixed composition at all energies. 
A simple linear fit, (X max ) = Diq ■ 
lg (E/eV) +c, yields an elongation rate of 54±2 
(stat.) gcm _2 /decade, but does not describe 
our data very well ( X 2 /Ndf= 24/13, P<3%). 
Allowing for a break in the elongation rate at 
an energy leads to a satisfactory fit with 
X 2 /Ndf= 9/11, P=63% and D 10 = 71 ± 5 
(stat.) gcm" 2 /decade below E h = 10 18 35 eV 
and Dio = 40 ±4 (stat.) gcm~ 2 /decade above 
this energy. This fit is indicated as a dashed 
gray line in Fig. 3. 
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Figure 4: Comparison of the measured elonga- 
tion rate, D\o, below (solid circle) and above 
(open circle) 10 18 35 eV to predictions of air 
shower simulations (red: protons, blue: iron). 

Due to the uncertainties of hadronic interac- 
tion at highest energies, the interpretation of 
these elongation rates is, however, ambiguous 
(cf. Fig. 4). Using the QGSJETII elongation 
rates the data suggests a moderate lightening 
of the primary cosmic at low energies and an 
almost constant composition at high energies, 
whereas the EPOS elongation rate is clearly 
larger than the measured one at high energies, 
which would indicate a transition from light 
to heavy elements. Theses ambiguities will be 
partially resolved by the analysis of the X max 
fluctuations as an additional mass sensitive pa- 
rameter. 

A comparison with previous measurements [9] 
is presented in Fig. 5. The results of all three 
experiments are compatible within their sys- 
tematic uncertainties. It is worthwhile noting 
that although the data presented here have 
been collected during the construction of the 
Pierre Auger Observatory, their statistical pre- 
cision already exceed that of preceeding exper- 
iments. 
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